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Abstract: The working group discussed several aspects of triple gauge coupling 
analysis viewed in the light of experiences with the first high energy data recorded 
at energies above the W pair threshold. Some analysis methods were reviewed 
briefly, and consideration given to better ways of characterising the data. The 
measurement of CP violating parameters was discussed. Results were prepared 
to further quantify the precision attainable on anaomalous couplings in the four- 
quark channel using jet-charge methods, and finally the trade off between maxi- 
mum LEP energy-vs-luminosity was quantified. 
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1 Introduction 



The triple-gauge-coupling (TGC) working group was composed of members of 
the theoretical and experimental community who are actively involved with the 
measurement and interpretation of TGCs at LEP2. Many of the participants also 
contributed to the LEP2 workshop in 1995. 

A review of the current status of the first experimental results from LEP2 was 
given in the plenary review [R.Sekulin - these proceedings] which also gives an 
overview of some of the calculational tools presently available as well as methods 
of analysis which have been used. 

At the LEP2 workshop many prophecies were made in advance of the first data 
taking above the W pair threshold. The group decided that the best use of this 
workshop would be to review how some of the ideas envisaged at that time had 
actually evolved in the light of experience with the first data, and to discuss 
some of the problems which have been brought to light. We also identified some 
quantatitive studies which were considered useful to present for future reference. 

In section two we consider briefly some of the limitations of present analysis 
methods. We then consider in some detail how experimental measurements might 
be better characterised in order to avoid, or minimise, some of these limitations 
in an unbinned maximum likelihood fit. Finally the possibilities for an improved 
fitting tool are described. 

In the third section we describe a development of a method employing optimal 
observables, which makes use of fully simulated distributions to extract TGC 
parameters. 

In the fourth section a study is presented in which the precision which can be 
obtained using the qqqq channel is quantified in the light of current experience 
of jet pairing and charge assignment purity. 

In section five we present a parametrization of the CP violating trilinear gauge 
couplings restricting ourselves to the gauge invariant operators with the lowest 
dimension. We also present some estimate of the sensitivity expected at LEP2. 

Finally in section six we present a study to quantify the variation of TGC mea- 
surement precision which would be obtained for different scenarios of maximum 
LEP2 energy versus the luminosity increase which results from sacrificing a few 
GeV from the absolute maximum energy which LEP can achieve. 

Throughout the following we make use of abbreviations for different W pair decay 
types. These are: (i) qqqq events where both W bosons decay to a quark- 
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antiquark final state, (ii) qq£z7£ events where one W decays to a quark-antiquark 
final state, and the other W decays to an electron, muon or tau, plus a neutrino 
and (iii) ^(Ivg events where both W bosons decay leptonically. We also refer to 
two distinct ways in which the phase space variables are normally characterised for 
each event. The first is to use quantities which are very close to those measured 
by the detectors, such as the four-vectors of the observed "fermions" (leptons 
and jets), this is referred to as 4-vector analysis. The second is to instead use the 
set of five angles given by the W~ flight direction with respect to the incoming 
beam direction, and the polar and azimuthal angles of the decay products of 
each W measured in the rest frame of the respective W. A full description of 
these angles can be found in [||, ^ These angles are commonly written as 
{cos 6*^^, cos 0*^, cos 6^*2, 0*2} and methods using these are referred to as 5- 
angle analyses. We refer to both schemes generically as the set {$}. 



2 Data characterisation for maximum likelihood 
fitting. 

Previous studies [|l| have examined possible observable effects from anomalous 
TGCs. Anomalous TGCs can affect both the total production cross-section and 
the shape of the differential cross-section as a function of the W~ production 
angle. Additionally, the relative contributions of each helicity state of the W 
bosons are changed, which in turn affects the distributions of their decay products. 

The most straightforward approach which has been considered is to use an un- 
binned multi-dimensional maximum likelihood method (ML) where one identifies 
the most likely measurement of a TGC a, by varying a to minimise the quantity 

-logL = - logP{^i,a) 

events i 

where $j represents the phase space measurements for event i. To do this is 
is necessary to have some function, P, which returns the probability density 
function for a particular $j to occur, and this can be conventionally derived as 
the normalised differential cross section. It is reasonably straightforward to code 
an analytic form for P if one restricts to on-shell formula without initial state 
radiation (ISR). The relevant expressions can be found in and an example 
code is given in |Q. Such approaches are simple, not very CPU intensive, but as 
well as neglecting ISR and Pw also neglect detector acceptance and resolution. 
The neglect of such effects was shown [|[| to lead to biases of roughly the same 
magnitude as the expected eventual experimental precision for 500pb~^. 
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More sophisticated approaches have been implemented, which attempt to allow 
for all of the above effects. A typical solution is to compare an observed data 
distribution with that predicted using fully simulated Monte Carlo events passed 
through the same selection requirements as the data. By using Monte Carlo 
distributions corresponding to different generated values of a one can make fits 
in order to measure the most likely true value of a. Examples of this are (i) a 
coarse binned maximum likelihood analysis [^] (ii) use of distributions of opti- 
mised observables p, (iii) fits to a reduced number of variables and (iii) 
nearest neighbour counting in Monte Carlo distributions [§]. These methods are 
straightforward in the case of one variable, but tend to require either a very large 
number of Monte Carlo events as the number of variables increases, or to ne- 
cessitate very coarse binning (the possible exceptions to this are the optimised 
observables methods described below). There are various ways to try to make 
these methods more efficient, however we decided to turn our attention back 
to the more straightforward maximum likelihood approach and consider how its 
failings could be rectified. 

In order to perform a "complete" ML analysis one would like to include all of the 
effects of ISR, Fw, detector resolution and detector acceptance into a likelihood 
function which we will call Pjuu- The issues involved in doing this fall into three 
classes: (i) to include ISR requires the inclusion of some internal integration over a 
radiated photon spectrum (ii) to include the effect of Fw depends upon whether 
one attempts to supply W mass information with each data event or instead 
wishes to internally integrate over the W virtualities and (iii) the facilities which 
can be included in Pfuii depend upon the choice of quantities which one measures 
for each data event. 

To completely specify an event requires 8 pieces of information. Each of the 
four decay fermions requires 3 quantities (masses are fixed). The total energy 
and momentum provide 4 constraints (assuming that the 4-momentum of any 
ISR photon is known). The total number of pieces of information required is 
therefore 4x3 — 4 = 8 (although the overall azimuthal angle of the event carries 
no information). In the 5-angle formulation the missing information is equivalent 
to ignoring the two W masses. It is also necessary to know the total energy and 
momentum to calculate the decay angles, and it is normally assumed that there 
is no ISR. This makes it impossible to then carry out a consistent integration over 
ISR in PjuU- Conversely if the 4- vectors of all of the decay fermions are specified 
then these completely determine both the W masses and the momentum carried 
by ISR. There is therefore no flexibility left for these to be integrated over in 

Pfull- 

Another important consideration is that if one chooses to characterise each event 
into 5-angles then all the quantities are complicated transforms of the basic de- 
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tector measurements such as lepton energy and direction etc. As a result the 
5-angles are highly correlated and rarely have Gaussian resolutions. Both of 
these facts make it non-trivial to incorporate integrations over detector resolu- 
tion or acceptance. If instead one chooses to characterise each event using the 
4- vectors then some of these limitations are removed. 

The group therefore considered the most desirable way to characterise each type 
of event from an experimental point of view. This was done separately for each 
W pair decay channel, and the results are summarised in table |l[ 

2.1 qqel^c and qq/il^^ events 

In these events the lepton is normally well measured and can be characterised 
by its energy, Ei, and its polar and azimuthal angles, cos 6*/ and (pi. These three 
quantities are largely uncorrelated and often have resolutions which may be ap- 
proximated [] as being Gaussian. 

The hadronically decaying W results in two jets. It is often (although not invari- 
ably) the case that the jet angle measurements are better than energy resolutions. 
It is therefore highly likely that any apparent jet mass is mainly determined by 
energy resolution rather than any physical mass. It therefore also follows that the 
measured invariant mass of the two jet system is poorly related to the true W vir- 
tuality on an event by event basis. This last point is important, as several fitting 
tools which exist take as input the pairs of measured 4-vectors belonging to each 
W, and then interpret the invariant mass of each pair literally when calculating 
the matrix element for the decay. This means that the inclusion of the effect 
of Fw is not via any integration internal to the fitting tool, but comes from the 
actual mass distribution of the data being fitted, and this is determined almost 
entirely by resolution errors. Although at this time no quantative study has been 
done, it seems unlikely that this is sensible. It is therefore desirable to have the 
option to allow Pfuii some freedom to integrate internally over W virtuality. We 
therefore concluded that the two jets should be characterised by their directions, 
cos^i, 01, cos 6^2 and 02 and optionally the sum of both jet energies, Ej^^^. 

This choice of these variables is designed to allow those that are well measured to 
be used, and those that may not we well measured to be integrated over easily. 
If -E/iad is used then 8 quantities are specified. Since it is not required to assume 
the total energy and momentum of the event to specify the 8 quantities, then the 
freedom is left for Pfuii to either assume no ISR or to integrate over a photon 

^ The energy will either be measured in a calorimeter or by the magnetic tracking detector, 
in which case either the measurement or its inverse is approximately Gaussian. 
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radiation spectrum. However there is no freedom left to integrate over the W 
virtuahties. If, however, it is considered that the hadronic energy measurement 
is relatively poor then one can either perform an explicit integration over the as- 
sumed resolution (see section |2.5|) or allow a complete integration over all possible 
values of Ej^^^ to be included in Pfuii- 



2.2 qqri'r events 

These events differ from the previous category because one observes only the tau 
decay products. It is likely that the direction of the decay products can serve 
as a useful approximation to the actual tau direction (with a larger resolution) 
but the energy will always be low as there is always at least one missing neutrino. 
Therefore the simplest way to specify the event is to use exactly the same quan- 
tities as for the qqeFg and qq/il^^ events, but in addition allow the integration 
over El in Pfuu- 



2.3 iuiiiyi events 

Studies have indicated that in the case of electron and muon events there is 
useful information on anomalous couplings. It is in principle possible to fully 
reconstruct [Q] these events up to a two fold ambiguity using the two observed 
leptons, the beam energy and momentum, and the W mass. However there are 
several problems in doing this due to experimental resolution and that one has no 
knowledge of the actual W masses for each event. The working group therefore 
considered that it would be better to use only the observed quantities and inte- 
grate over everything else in Pfuii- Thus such events would be characterised by 
Ei+ , cos 9i+ , (f)i+ , El- , cos 9i- and 0^- (note again that only the relative carries 
information). If in addition an integration over Ei+ and Ei+ can be included then 
tau events could also be used. 



2.4 qqqq events 

These are the most problematic. The experience of the working group from 
various Monte Carlo studies is that the W production direction can be well mea- 
sured from these events, but as for the qqtp£ events there is probably no useful 
information about the individual W masses on an event by event basis. However, 

^ It was not considered sensible (at this stage at least) to consider incorporating an integra- 
tion over all tau decay modes into Pfuii ■ 
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following the previous reasoning, we might think of specifying all of the jet angles, 
cos^j, 0i, i = 1,4:. In this case we would have given all 8 pieces of information, 
hence completely specifying the W masses, once the total energy and momen- 
tum are fixed (in fact this is used by DELPHI in one of their mass measurement 
methods) . At the present time no detailed studies have done to resolve this point 
from the point of view of TGC analysis, and further work is required. 

2.5 Inclusion of detector resolution and acceptance 

The quantities chosen to specify each type of event are largely uncorrelated and 
are often described by Gaussian resolutions. As such it is possible to imagine 
performing a simple Gaussian integration over each to allow for detector resolu- 
tion. This can be done externally to Pfuii using commonly available integrating 
functions. 

Exactly how many quantities to integrate over depends upon the event type. In 
the case of o^iPe events it may be that the precision of the lepton measurements is 
more than adequate for no resolution on either energy or direction to be included. 
In this case one needs to perform an integration over only the jet angles. By 
the same arguments it may be that in the case of tUglui events resolution can 
be neglected entirely. Exactly how to proceed in each case is an experimental 
question, and also depends upon the statistical precision of the data. 

Similar approximations may be used in order to account for detector acceptance. 
For instance detectors normally have a sharp cut off in cos^ for lepton iden- 
tification, but otherwise are near 100% efficiency and uniform in 0. A simple 
approximation such as this can be included in the normalisation of Pfuii- In the 
case of jets the situation is not always so simple, but it may nevertheless be ap- 
propriate to apply a sharp fiducial cut such that jets are only used in a region 
where they are well reconstructed with high efficiency. 

The treatment of resolution and acceptance suggested here is an experimentally 
motivated approximation, and as such will never be "exact". However since 
the biases incurred by completely ignoring these eff'ects are only expected to be 
approximately the size of the final LEP2 statistical error, it may be that the level 
of approximation involved here is adequate for all practical purposes. Further 
studies are needed to quantify this issue. 
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2.6 Prospects for a new semi-analytic analysis tool 

A feasibility study of a semi-analytic analysis tool has been made. For this pur- 
pose two different problems should be dealt with. In the first place an efficient 
program should be developed which for each set of variables from Table 1 cal- 
culates the probabilities. In that evaluation the effects of the W-width and ISR 
should be taken into account. Secondly, it should be studied whether this an- 
alyzer properly determines the TGC parameters from a Monte Carlo sample of 
unweighted events. 

The first issue, developing a program which calculates cross sections differential 
in the appropriate variables has been accomplished. In fact, two independent 
programs have been made in order to have cross checks. One program uses the 
matrix elements of ERATO, the other of EXCALIBUR. All kinematical variables 
of table 1 have been implemented, except the case with the tau-lepton. Without 
ISR the chosen variables require two, one integration, sometimes none. With ISR 
two more integrations have to be carried out. 

The second issue, whether maximum likelihood fits based on probabilities related 
to the above specific kinematic variables can really determine TGC parameters 
is at present under study. One generates unweighted events, uses of these events 
only the variables of table 1 and calculates the likelihood function for e.g. nine 
values of a specific TGC parameter. From this the TGC parameter is determined. 
In this way one could try to determine any parameter on which the matrix element 
depends. 

In fact, so far the most extensive studies have been made for a parameter other 
than a TGC. The W-mass has been chosen as first parameter to be studied. One 
reason for this choice is that the question of experimentally motivated quanti- 
ties can also be relevant for the mass determination. Another is, that one may 
eventually like to determine simultaneously both the mass and a TGC coupling. 

So, as a first test of the method the mass determination from experimentally 
motivated quantities was chosen. To this end typically 1600 unweighted events 
were generated with ERATO or EXCALIBUR. The semi-leptonic case has been 
studied in great detail. The data samples are CC03, CCIO and CC20 based sets 
including ISR and have been analyzed with the same matrix elements or simpler 
ones (e.g. a data set with ISR analyzed without ISR). 

The following results can be reported. Within the errors, typically 35-50 MeV, 
the input mass could be reconstructed when the same matrix element is used in 
the analysis as in the generation of the sample. When simpler matrix elements 
in the analysis than in the generation are used, shifts can arise. Large shifts (300 
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MeV in some cases) may occur due to the neglect of ISR in the analysis. Smaller 
shifts (up to 50 MeV) are found when a CC20 sample is analyzed with CC03 
probabilities. 

The mass determination is also successful in the four quark case when one assumes 
to know which quark pair comes from a W. This reduces the folding. When the 
method is applied to the purely leptonic case the mass reconstruction is not 
efficient, as one would expect as there is less kinematical information. 

The results on the mass determination will be published in more detail |^ else- 
where. It seems that the maximum likelihood method to determine Mw from 
experimental motivated quantities is a useful addition to the existing direct re- 
construction method. 

As to the TGC determination, the first results look promising, but more results 
should be obtained before a detailed discussion of the merits of the method can 
be presented. However it already now is clear that the method takes ISR and full 
matrix elements into account in a very satisfactory way. 



Channel 


Well measured 


"Poorly measured" 
(may be integrated over) 


qqeI7e 


El, cos 01, (pi, cos 6i, 01, cos 62, 02 


rptot 
^had 


qq^z7^ 


// 


II 


qqTz7^ 


COsOl, (f)i, COS 6^1, 01, COS 6^2, 02 


rp rptot 
^had 




Ei+ , COS 9i+ , 0;+ , El- , COS 9i- , (f)i- 




qqqq 


cos6i,(f)i, i = l,4 





Table 1: Experimentally motivated quantities to specify for different W pair 
event categories. 



3 Development of optimal observable analysis. 

In section |^ several methods of TGC analysis were mentioned which employ fully 
simulated Monte Carlo events in order to analyse the data. These methods 
include all effects of ISR, Fw, detector acceptance and resolution. One problem 
which often arises is the number of bins necessary in a multidimensional analysis. 
However methods which use "optimal observables" have the potential to avoid 
this problem. 

The use of optimal observables (00) has been discussed widely in the literature [|I], 
pUj , pT] , . The method consists in constructing observables, namely functions of 
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the observed kinematic variables $, which maximise the sensitivity to anomalous 
couplings There is one observable per coupling, Cq,($), given by: 



= £ (1) 

where the functions ^ are obtained from the quadratic form of the Born level 
cross section: 

da dao dai nda2 , 

— = ha ha (2) 

d^ d^ d^ d^ ^ ' 

The observable is calculated for each event from the measured $ variables. In 
principle the mean value of O (averaged over the data set) contains all of the 
information which can be extracted from the $ distribution of the data, provided 
that it is small enough that the neglect of the second order term has a negligible 
effect. The spectrum of the optimal observables might be regarded as a projec- 
tion of the five- dimensional differential cross section in $ onto a direction which 
optimises the sensitivity. 

As originally envisaged and discussed in the literature, the mean value is used 
to determine a through the inversion of a matrix relation between the two. The 
effects of ISR and Fw can only be included if the the resulting functional form 
of the differential cross section of equation |I| is known, which is not generally the 
case, and therefore biases are introduced. It is also in principle possible to include 
detector effects in the procedure, but this can lead to problems of instability. 

Instead of following this procedure a recent OPAL analysis adopts a different 
approach. A similar method has also been developed by DELPHI |p!3|. 



In the OPAL analysis each event is analysed using the 5-angles as the phase 
space variables. Each 00 quantity is calculated using the the on-shell/no ISR 
differential cross section (it is straightforward to obtain an analytic function for 
this). The value of the 00 measured per event is then treated just as any 
other experimental quantity whose distribution is sensitive to variations in the 
relevant anomalous coupling. As example the observed distribution of the optimal 
observable for the a\/i/^ parameter is shown in figure |I] for data and Monte Carlo. 
The expected distribution for all values of a is obtained using fully simulated 
EXCALIBUR events generated for a few specific values, and then using the 
quadratic dependence upon a in each bin of the spectrum of Cq,. The Monte Carlo 
distributions are generated with ISR and Fw switched on. The contribution from 
background contamination is added to the spectrum. 

To measure the TGC parameter a binned likelihood fit of the observed data to 
the expected distribution predicted by the Monte Carlo is performed. Since like 



10 



is compared with like the method includes all of the required effects and no bias 
should be introduced. 

There are two drawbacks of the method. The first is that the 00 is calculated 
using an on-shell/no ISR differential cross section. This can be considered as an 
approximation to the true 00. The approximation does not introduce any bias, 
nor should it be confused with the fact that the data is fitted to Monte Carlo 
which includes ISR and Fw, however it does in principle lead to loss of optimality 
which means that the measured statistical error will increase. It is for this reason 
that the full distribution of the calculated for each event is used, making 
better use of the available information than the mean value in the case of this 
non-optimality. The second drawback is that the simple 00 method is based 
upon an expansion about a single value of a given TGC and is only optimal if 
the fitted values lie in a small range around this value. At present the statistical 
precision of the data is rather poor and so the 00 method will not necessarily give 
the true confidence levels for large values of the errors. However this limitation 
is expected to diminish as the statistical precision of the data increases in future 
years. However this method is in practice almost adequate to analyse the current 
small data set. The logL distributions obtained from the method described, and 
from a more conventional binned maximum likelihood analysis of the same data 
give almost identical errors and central values. These log L curves can be seen in 
reference 0. 



4 The use of jet charge in qqqq events 

In events where both W pairs decay hadronically the simple method of deter- 
mining the W charge used in the semi-hadronic and leptonic events is no longer 
available. A method of assigning a charge to the Ws from the properties of the 
jet pairs into which they decay must be sought. Various possibilities suggest 
themselves, but none of them lead to identification of the W charges with 100% 
certainty. Previous studies have either assumed that the charge is unknown or 
completely determined; here we investigate the effect that varying the efficiency 
of the W charge determination has on the ability to measure aw<f,- We generated 
events with Pythia. The correct pairing of jets into Ws was taken and, for a 
fraction (l — fg) of the events (chosen randomly), the correct charge assignment of 
the Ws was taken. For the remaining fraction fg of the events, the pair constitut- 
ing the W~ was called W+ and vice-versa. This procedure defines the "assigned 

^ Alternately, an extension of the method is possible whereby the term is incuded. 

* If the measured central value, ao, is different from the expansion point of a = then the 
method can be iterated to use 00s calculated from a functional form of equation ^ where a is 
replaced by a — ag- 
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Figure 1: The distribution of the optimal observable for awij> obtained from the 
qq£z7£ events. The hatched histogram shows the non WW background. These 
are compared with the distribution expected in the Standard Model using fully 
simulated Monte Carlo events. The predicted distributions for aw(i> = +2(— 2) 
are also shown as dotted (dashed) lines. 

momenta" with a probability fg of assumed incorrect jet charge assignment. 

In the analysis of the events, the unnormalized probability for observing each 
event was taken as (1 — /a)p(l,2) + /a(2, 1), where p{i,j) is a squared matrix 
element. For = (1, 2) the production and decay angles of W~ and W"*" were 
calculated using the "assigned momenta" defined above, while for = (2, 1) 
the production and decay angles were taken with the "assigned" W~ interpreted 
as the W"*" and vice-versa. 

Results of extended maximum likelihood fits of aw(f> to 1000 events at 190 GeV 
are shown in table |^. The decay angles from both Ws were included in the fits, 
but folded, thus assuming no quark flavour information in the W decay products. 
Thus each p{i,j) is an average of four matrix element calculations. 

There is a definite pattern in the results: 

• If the correct misassignment probability is assumed in analysis (the diagonal 
entries in the table), the result is unbiased, but the precision decreases 
rapidly with increasing fg. For fg = 0.5, the result is identical with that 
obtained simply by folding production and both pairs of decay angles in the 
fit (an eightfold ambiguity), i.e. without applying any jet charge selection: 
(Xw^ = 0.024 ± 0.064. For fg = 0.2 (a value which might be expected in 
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fa, Rate 
assumed 

in analysis 


fg, Misassignment rate applied to generated events 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.00 


0.004 
±0.039 


-0.189 
±0.033 


-0.298 
±0.030 








0.05 




-0.095 
±0.040 










0.10 




-0.011 
±0.048 


-0.129 
±0.042 








0.15 




0.057 
±0.055 


-0.053 
±0.049 








0.20 




0.103 
±0.058 


0.019 
±0.055 


-0.081 
±0.051 






0.25 






0.056 

±0.051 








0.30 






0.084 
±0.061 


0.0255 
±0.060 


-0.051 
±0.057 




0.40 








0.060 
±0.063 


0.023 
±0.063 


-0.020 
±0.061 


0.50 










0.024 
±0.064 


0.024 
±0.064 



Table 2: Results of fits to aw4, 
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practice), the precision has already deteriorated to a value rather close to 
this. 

• If an incorrect misassignment probability is assumed in analysis, the result is 
biased. For instance, if the true misassignment rate is fg — 0.2, but a value 

fa = 0.0 is taken in analysis, a large negative bias in a^/^f, is induced. This 
arises from the fact that the jet charge misassignment has had the effect of 
flattening the production angular distribution, which is also predicted for 
negative values of aw4>- Unfortunately, in the region fg = 0.2, the bias seen 
in the table varies rather rapidly as the assumed misassignment probability 
diverges from its true value, implying that a large systematic error might 
be incurred in the fitted TGC value. 

These considerations will have to be taken into account in determining the use- 
fulness of jet charge information in the analysis of the qqqq channel. 

This above analysis assumes that the angle of the W is well measured. We then 
wished to answer the question "What effect will imperfections in our definition 
of the W direction have on our ability to measure aw<i) ?" 

The standard of way of estimating the W direction is to divide the tracks (and 

energy clusters) of the event into 4 "jets", and to add the momenta of pairs of 
jets together to make the W momenta. The resulting directions are not aligned 
precisely with the directions of the Ws. 

In events where the reconstructed W is badly misaligned, the discrepancy can 
sometimes be traced to a misallocation of the jets to the Ws. This leads naturally 
to the idea of a measurement degradation due to a jet misassignment. However, 
this is an oversimplification. In another class of badly measured events it is the 
"jets" themselves which fail to measure the direction of the fundamental decay 
fermions. To be more precise, in the context of Pythia, if the sets of tracks 
which make up the four jets is compared with the sets of particles which are the 
decay products of the four fermions, then for these events one or more of the 
reconstructed jets will not be made up of particles from the decay of a single 
fermion. 

The details of the problems caused by these effects may depend rather sensitively 
on the apparatus used in the measurement and the algorithm used to extract 
the W directions. For this reason we made no attempt at the sort of generic 
measurement which is suitable for the effects of jet charge. Instead, in order to 
provide an idea of the sort of degradation which follows from these considerations, 
we have taken the example of the OPAL apparatus and processed EXCALIBUR 
events through a a full simulation of the detector. The resulting simulated data 
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is divided into 4 jets using a Durham type algorithm; a kinematic fit is performed 
on jets constraining the pairs to have the same mass, and the pairing with the 
highest mass is taken as the correct one. The charge of each jet is defined using a 
momentum weighted charge measure and the pair of jets which have the highest 
summed jet charge is taken as the positive W. For a dataset consisting of 7,000 
fully simulated qqqq events we fit the value of aw<t, for all the events which are 
successfully reconstructed as four jet events. This defines a 'realistic measure' 
of our ability to measure aw(t>- Then the opening angles flj^ between each jet 
and fermion are calculated, and the subset of events is taken where the minimum 
value of fljtj, allows a unique association between quark and jet. This looses about 
10% of the data sample. For these events the W's direction is calculated by taking 
the correct pairing as given by the generator information and from that a value 
of cos^u; is derived. The charge assignment is then made by the normal charge 
weighting methods. It turns out that for these events the alignment between jet 
and quark is good, but the jet energies are not nearly so tightly constrained. Thus 
the W direction is much less well measured than the individual quark directions. 
These are the "good jet matching events" . Finally we used the correct value of 
the W charge, and a value and error was calculated for these "charge correct" 
events. 



Data 
Set 


Reconstructed 


Pairing 

correct 


Charge 
correct 


aw4> 


0.03-0.077+0.082 


-0.18-0.072+0.080 


0.00+0.028 



Table 3: Results of fits to aw(p 



The values of aw<j> which were found can be seen in Table |^, where it can be seen 
that the dominant cause of measurement degradation is the mismeasurement of 
the W charge. Reanalysing the 90% sample without any generator information 
produces a value fully consistent with the original sample. 

These values should not be extrapolated to other values of aw^- The accuracy of 
the measurement of aw(j> depends on the value of aw(t>, and experimental effects 
add towards that total error in ways which may also vary with aw(j>- 



5 CP violating TGC parameters 



In the usual phenomenological parametrization of the TGC |]1| the CP violating 
interactions are also included in the form: 



^TGC — e Qvww 



15 



w 



(3) 



m 
where 

and 

In Eq.(|D is the ly-boson field, and the usual definitions g^ww = 1, dzww = 
ctg 9w are used. Finally e = g sm9w = g' cos9w 

What is the meaning of the above Lagrangian ? One way to answer this ques- 
tion is to look at it as the low energy limit of a manifestely SU{2)l x U{1)y 



gauge invariant theory. This can be done |T^, [T^, in an effective Lagrangian 
approach by considering gauge-invariant operators involving higher-dimensional 
interactions among the gauge bosons and the Higgs field. These operators are 
scaled by an unknown parameter K^p describing the characteristic scale of some 
high energy New Physics, generating at low energies the effective interaction Ltgc 
as a residual effect. In order to generate all kinds of TGC appearing in Eq.(|D, 
we need operators with dimension up to twelve. On the other hand, restricting 
ourselves to SU{2)l x ?7(l)y- invariant operators with dimension six, which are 
the lowest order ones in a l/A^vp expansion, we end up with the following list of 



operators capable of inducing CP-violating TGC couplings [IS 



Ow = l^iW^p X W^^) ■ W\ , (4) 

where 

B'^'' = ^E'''"'^Bp^ , ly'"' = le/^-P-iy^, (5) 
and $ is the Higgs doublet. 

The New Physics contribution from the above operators is described by the ef- 
fective Lagrangian 

gg' asw . c^w /^x 

where the relations between asw^ ctw, and the parameters appearing in the 
Eq.(|^) are given by 

kz = — tan^ Oyjk^ \z = X-y ■ (7) 



16 



In order to get an estimate of the expected sensitivity at LEP2, we performed 
an analysis based on the event generator ERATO, using the method of Optimal 
Observables, but without taking into account any experimental effect such as 
resolution and related background. The results are given in Table § corresponding 
to an integrated luminosity of 500 pb~^. 



(GeV) 


161 


172 


192 


205 


OtBW 


1.81 


0.79 


0.08 


0.08 




0.43 


0.14 


0.02 


0.016 



Table 4: One standard deviation errors on TGC parameters. 

It should be emphasised that the limits on the CP violating couplings are expected 
to be of the same order as those on the CP conserving ones. 



6 High energy -vs- high luminosity 

The maximum centre-of-mass energy attainable at LEP is expected to be limited 
by the radio frequency accelerating gradient available. As a result, a trade-off 
will be inevitable between the maximum centre-of-mass energy achievable and the 
reliability of machine operation, by keeping a fraction of the RF in reserve. Other 
machine performance issues emphasise this point further. There are therefore 
a range of strategies for running LEP during the LEP2 programme. The two 
extreme cases may be characterized as: 

1. run at the highest energy attainable, with reduced safety factors, and accept 
the loss in integrated luminosity; 

2. run at a lower energy and attempt to collect as much integrated luminosity 
as possible. 

In order to seek a balance between these extremes, one should take into account 
the extra reach for searches and the expected precision of measurements. One 
element in this accounting is the variation of the expected precision for TGC 
measurements as a function of the centre-of-mass energy and the integrated lu- 
minosity. 
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6.1 Sensitivity studies with simulated events 



The prospective sensitivity to the anomalous couphngs has been studied for the 
case of the a parameters in the W^, and W models. Two independent 
event generation and analysis procedures were applied and consistent results were 
obtained. 



In the first analysis, events were generated with the EXCALIBUR ||T4| program 
at various a values, and analysed with a fitting program using the ERATO |]19 
matrix element routines to calculate the differential cross-section. The second 
analysis employed the WOPPER pO|] event generator, and analysed events using 
the on-shell cross-section formulae from Bilenky et al. [§. These two analysis 
strategies differ in sophistication, and include slightly different information about 
the anomalous coupling sensitivity but for the purpose of this study, i.e. the 
determination of the sensitivity to the anomalous couplings, they give similar 
results. 



Following [0, only the semileptonic channels qqei^e and qq/iz/^ were considered, 
and no detector simulation was applied. These "ideal detector" studies do not, of 
course, correspond directly to real analyses of events observed in the LEP detec- 
tors, because various additional effects enter: limited acceptances, backgrounds, 
and detector resolution degrade the performance. On the other hand, including 
information from the qqrz/,-, \^v\~v and qqqq channels improves the precision 
relative to the two-channel ideal detector case. 

It is possible to estimate the statistical precision that might be obtained in 
real analyses compared with these idealized studies by comparing the precision 
achieved and expected from the small 172 GeV data sample already taken. The 
calibration factors were found to be within 20% of unity within most cases, sup- 
porting the usefulness of the idealized studies. Only statistical errors were consid- 
ered in these studies. While systematic effects may eventually become significant 
for these measurements, the dominant systematic contributions are expected to 
come from detector effects, which can presumably be understood with sufficient 
precision with sufficiently large data samples. 

Figure |^ shows the predicted error on the a parameters for centre-of-mass energies 
ranging from 170 GeV to 210 GeV, for 500 pb~^ collected by a single experiment. 
The errors are derived for the case where the central value of the a parameter is 
close to zero — it was found that the expected error decreases by (9(30%) if the 
true a parameter is 1, for example, but the main effect is simply the substantially 
increased W-pair production cross-section for such non-standard couplings. The 
sensitivities at four canonical centre-of-mass energies are shown in table |^. It 
is evident that the gain in precision with energy is very marked below around 
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180 GeV, but that the gain then slows, although the size of this effect differs 
for the different couplings. It is interesting in particular to compare the 190 and 
200 GeV points: a unit of luminosity at 200 GeV is approximately as sensitive 
as 1.4 to 1.8 units at 190 GeV. 



6.2 Conclusion 



The best measurement of anomalous coupling parameters at LEP2 requires that 
a large fraction of the luminosity should be taken above 180 GcV centre-of-mass 
energy, as is planned. It is important that the luminosity available should be 
high: the 500 pb^^ considered in previous studies remains an apposite goal. The 
sensitivity per unit luminosity rises above 180 GeV, so it is useful to increase the 
ccntre-of-mass energy if only a modest cut in luminosity is taken. The sensitivity 
rise is, however, slow, so this analysis disfavours operating strategies where a 
substantial luminosity penalty is incurred for the sake of only a few GeV in 
centre-of-mass energy. 



7 Summary 



The first data which has been recorded at LEP2 above the W pair threshold has 
been followed by a rapid analysis by all of the LEP experiments to search for 
possible signals for anomalous triple gauge boson couplings. That this could be 
done so quickly is a tribute the the large amount of preparatory work performed 
by theoreticians and experimentalists in the provision of tools, and the develop- 
ment of analysis methods. In this working group we reviewed our experience in 
the light of this first data. We primarily discussed the limitations which arise 
when trying to incorporate all effects of W width, ISR, detector acceptance and 
resolution into an analysis. These effects had been previously shown to lead to 
biases, and several techniques have been developed to try to handle this. Most 
methods do however have some limitations, and so we have considered alternative 
ways to characterise the data which may be more amenable to the use of analytic 
fitting tools. One of the most important outcomes of the work was to nucleate 
the provision of a new fitting tool. Work on this is still underway and will be 
published separately at a later date. Another important topic is the extent to 
which use can be made of W pair events decaying into the four quark channel. 
Analysis of this channel is difficult due to the problems in identifying the correct 
jet pairing and the correct W charge assignment, and therefore a study was made 
to quantify the degradation in precision due to these effects. We also briefly 
discussed the formulation of CP violating TGC parameters in a gauge invariant 
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Centre-of-mass energy 




AaBct> 




172 GeV 


0.042 


0.25 


0.067 


183 GeV 


0.027 


0.12 


0.042 


190 GeV 


0.023 


0.092 


0.035 


200 GeV 


0.020 


0.068 


0.029 



Table 5: Expectations for triple-gauge coupling precisions from a single experi- 
ment, approximated by the "ideal" precision for the qqez/g and qq/if^ channels. 
In each case the precision which would be obtained with 500 pb^^ at that centre- 
of-mass energy is given. 
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Figure 2: Evolution of the expected statistical precision on the measurements 
of the anomalous coupling parameters a for different models, as a function of 
centre-of-mass energy. The precisions given correspond to a single ideal detector 
analysis of the qqei^e and g^ht^h channels, employing the production and decay 
angles of both W's, but without resolving the two-fold decay ambiguity from the 
hadronically decaying W. 
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framework. Finally we have presented a study of the trade off between achieving 
the highest possible LEP2 energy at the expense of integrated luminosity. 
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